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The interplay between disorder, quantum fluctuations and dissipation is studied in the random transverse Ising 
chain coupled to a dissipative Ohmic bath with a real space renormalization group. A typically very large length 
scale, L* , is identified above which the physics of frozen clusters dominates. Below L* a strong disorder fixed 
point determines scaling at a pseudo-critical point. In a Griffiths-McCoy region frozen clusters produce already 
a finite magnetization resulting in a classical low temperature behavior of the susceptibility and specific heat. 
These override the confluent singularities that are characterized by a continuously varying exponent z and are 
visible above a temperature T* ~ L*^^. 
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The presence of quenched disorder in a quantum mechani- 
cal system may have drastic effects in particular close to and at 
a quantum critical point. The appearance of Griffiths-McCoy 
singularities leading to the divergence of various quan- 

tities like the susceptibility at zero temperature even far away 
from a quantum criticalpoint, has received considerable atten- 
tion recently r3l |4l yl la, Iztl ■ This quantum Griffiths behavior 
is characteristic for quantum phase transitions described by an 
infinite randomness fixed point (IRFP) j^, which was shown 
to be relevant for many disordered quantum systems |^. 

Quantum Griffiths behavior was proposed to be the physical 
mechanism responsible for the "non-Fermi-liquid" behavior 
observed in many heavy-fermion materials llOtll ill . However, 
it was also argued that in a dissipative environment, as in met- 
als due to the conduction electrons, such a quantum Griffiths 
behavior might essentially be non-existent fl2l [Till . More- 
over, even the underlying sharp quantum phase transition it- 
self was shown to be rounded in dissipative model systems 
Obviously there is a need to examine carefully the ef- 
fect of dissipation on a quantum system displaying IRFP and 
quantum Griffiths behavior in the non-dissipative case and to 
treat correctly the mixing of critical and Griffiths-McCoy sin- 
gularities — which is what we intend to do in this letter. 

The properties of a single spin coupled to a dissipative bath 
has been extensively examined 1 15]. Upon increasing the cou- 
pling strength between spin and bath degrees of freedom it 
displays at zero temperature a transition from a non-localized 
phase, in which the spin can still tunnel, to a localized phase, 
in which tunneling ceases and the spin behaves classically. 
Such a transition is also present in an infinite ferromagnetic 
spin chain coupled to a dissipative bath, as it was recently 
shown numerically |16]. Here we want to focus on the in- 
terplay of disorder, quantum fluctuations and dissipation and 
study the Random Transverse Field Ising Chain (RTIC) where 
each spin is coupled to an ohmic bath of harmonic oscillators 
ivh . It is defined on a chain of length L with periodic bound- 
ary conditions (p.b.c.) and described by the Hamiltonian H: 



L 



Pi. 



(1) 



where (S]'~ are Pauli matrices and the masses of the oscilla- 
tors are set to one. The quenched random bonds 7/ (respec- 
tively random transverse field hi) are uniformly distributed 
between and Jq (respectively between and ho). The 
properties of the bath are specified by its spectral function 
-/<(«) = f UC|,/(0^,5(co-(0^,) - fa,coe-"/"' where i2, is 
a cutoff frequency. Initially the spin-bath couplings and cut- 
off frequencies are site-independent, i.e. a, — a and i2, = H, 
but both become site-dependent under renormalization. 

To characterize the ground state properties of this sys- 
tem Q, we follow the idea of a real space renormalization 
group (RG) procedure introduced in Ref. 1 18] and pushed fur- 
ther in the context of the RTIC without dissipation in Ref. 1 3]. 
The strategy is to find the largest coupling in the chain, either a 
transverse field or a bond, compute the ground state of the as- 
sociated part of the Hamiltonian and treat the remaining cou- 
plings in perturbation theory. The bath degrees of freedom are 
dealt with in the spirit of the "adiabatic renormalization" in- 
troduced in the context of the (single) spin-boson (SB) model 
fTill . where it describes accurately its critical behavior 1 1^. 

Suppose that the largest coupling in the chain is a trans- 
verse field, say /12. Before we treat the coupling of site 2 
to the rest of the system —Jia\a2 — J2O2OI perturbatively as 
in 1 3] we consider the effect of the part —h2<J^ + L<:(p^/2 + 
(Okxj_ /2 + CkXk<3\ ) of the Hamiltonian, which represents a sin- 
gle SB model. For this we integrate out frequencies co^. that 
are much larger than a lower cut-off frequency ph2 <C ^2 with 
the dimensionless parameter p':^ \. Since for those oscilla- 
tors co^ ^ h2 one can assume that they adjust instantaneously 
to the current value of ol tiie renormalized energy splitting is 
easily calculated within the adiabatic approximation lisll and 
one gets an effective transverse field h2 < h2: 



h2 = h2 (p/l2/n2) 



«2 



^2 — ph2 . 



(2) 
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If h2 is still the largest coupling in the chain the iteration (|2} 
is repeated. Two situations may occur depending on the value 
of a2- If a2 < 1 this procedure will converge to a finite 
value hi = h2{pli2/^2)"'^^''^^'^^ and the SB system at site 2 
is in a delocalized phase in which the spin and the bath can be 
considered as being decoupled (formally a2 = 0), as demon- 
strated by an RG treatment in 1 19|. 

If this value hj is still the largest coupling in the chain it 
will be aligned with the transverse field. As in the RTIC with- 
out dissipation, this spin is then decimated (as it will not con- 
tribute to the magnetic susceptibility) and gives rise, in second 
order degenerate perturbation theory, to an effective coupling 
Ji between the neighboring moments at site 1 and 3 1 3j| 



/l = JiJ2/h*2 



(3) 



If a2 > 1, h2 can be made arbitrarily small by repeating the 
procedure (|3 implying that the SB system on site 2 is in its 
localized phase [19] and essentially behaves classically: the 
decimation rule (|2j indeed amounts to set h2 = 0. Such a mo- 
ment, or cluster of spins, will be aligned with an infinitesimal 
external longitudinal field and is denoted as "frozen". 

Suppose now that the largest coupling in the chain is a 
bond, say 72- The part of the Hamiltonian that we focus on is 
-720503 +L/=2,3 Lit,- (/'^y2 + 0)A.,.xfy2 + Q,.XA-,.oj), i.e a sub- 
system of two spin-bosons coupled via J2- We find that in 
second order perturbation theory the ground state of this sub- 
system is equivalent to a single SB system coupled to both 
baths leading to the additive rule 



^2 = 0C2 + as 



(4) 



Integrating out the degrees of freedom of both baths with fre- 
quencies larger than pJ2 (as done previously for the single SB 
system |19]) the two moments at 2 and 3 are replaced by a 
single moment p(2 with an effective transverse field h2: 



P2 



' J2 {^2/ {^3 J 
:/J2+A'3 , Cl2^pJ2 



(6) 



where /j; is the magnetic moment of site / (in the original 
model, one has /j, = 1 independently of /)• Combining Eq. (|4} 
and Eq. (|5Jl one clearly sees that a,- ~ /J,a. 

In the following we analyze this RG procedure defined by 
the decimation rules ( I2I6> numerically. This is done by con- 
sidering a finite chain of size L with p.b.c. and iterating the 
decimation rules until only one site is left. This numerical im- 
plementation has been widely used in previous works |9], and 
it has been shown in particular to reproduce with good accu- 
racy the exact results of Ref. 1 3] for the RTIC. We fix ho = 1 
and concentrate on the parameter plane {a,Jo). All data were 
obtained by averaging over 10^ different disorder realizations 
(if not mentioned otherwise), and the disorder average of an 
observable O is denoted by 0. The decimation rules (I2I6> 
depend explicitly on the "ad hoc" parameter p (or more pre- 
cisely on the ratio Q./p). For the moment we fix = lO'* 
and discuss the weak dependence on this parameter below. 
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FIG. 1: a) PL(log(ro//i)) as a function of log(r()//j) for different 
system sizes L for a = 0.03 (/iq = 1, = 0.34), i.e. far from the 
pseudo-critical point, b) P/.(log(ro//t)) as a function of log (ro//iL^) 
for different system sizes for a = 0.03 (/jq = I, Jq = 0.34). The 
straight dashed line has slope l/zwithz= 1.65(5). c) ^L(log(ro//j)) 
as a function of log(ro/A) for different system sizes L for a = 
0.052 (h{) = I, Jo = 0.34), i.e. at the pseudo-critical point, d) 
L^Pz. (log (To//!)) as a function of log(ro//2)/L*f for different sys- 
tem sizes for a = 0.052, \\i = 0.32. 



The transverse field h acting on the last remaining spin is 
an estimate for the smallest excitation energy. Its distribution, 
Piih/ro), where To is the largest coupling in the initial chain 
of size L, reflects the characteristics of the gap distribution [6J. 
Since the last spin can either be frozen (i.e the last field h is 
zero) or non-frozen we split PL{h/YQ) into two parts: 



PL{h/U) ^AlPlQi/To) + (1 -Ai)5(/2/ro) 



(7) 



where Pi{h/Y()) is the restricted distribution of the last fields 
in the samples that are non-frozen and Ai is the fraction of 
these samples. It, or equivalently Pi(log(ro//i)), represents 
the distribution of the smallest excitation energy in the en- 
semble of non-localized spins. 

Let us first present data obtained for Jq — 0.34. Fig. 
shows Pz,(log(ro//!)), for a — 0.03. For a system close to, 
but not at, a quantum critical point described by an IRFP one 
expects indications of Griffiths-McCoy singularities charac- 
terized by the following scaling behavior for Pl 



PL{\og{Yalh)) = 'P{\og{YQ/hU)) , 



(8) 



where z is a dynamical exponent continuously varying with 
{Jq, a, etc.). Fig. [Qi shows a good data collapse with z = 1 .65 
for the chosen coupling constant a = 0.03. The slope of the 
dotted line in Fig. ^ is identical to 1 /z and upon increasing 
a we observe that the slope, 1 /z, decreases. Our numerical 
estimates for l/z(oc) are shown in Fig. |2j they indicate that 
l/z approaches zero for some critical value tte, which implies 
formally z °° for a — > Uc. This is also characteristic for an 
IRFP, where PL(log(ro//i)) is expected to scale as 



PL(iog(ro//i)) =L-^fPiRFp(L"^iog(ro//i) . 



(9) 
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FIG. 2: Magnetization mgq and inverse dynamical exponent 1/z as a 
function of a (for ho = 1, io = 0.34) 




FIG. 3: Magnetic moment a/j(L) as a function of the system size L 
for different values of 5. Inset : (X/7(L) for L > 10^ as a function of 
L and for A = 0.13. The linear behavior suggests a finite nieq. The 
lines are guides to the eyes. Due to the large system sizes the data 
are averaged only over 500 different disorder realizations.) 



\|/ is a critical exponent characterizing the IRFP. Fig. shows 
FL(log(ro//i)) for a — 0.052: one observes that it broadens 
systematically with increasing system size, in contrast to the 
data shown in Fig. and Fig. [Qj displays a good data 
collapse according to (|5Jl with \\t = 0.32. Varying a only 
slightly worsens the data collapse substantially, therefore we 
take = 0.052 as our estimate for the putative critical point 
(for ho = 1 and Jo = 0.34) and denote by A = (a^ - a) /a^ the 
distance from it. 

The magnetic moment /j of the last remaining spin in the 
decimation procedure represents an estimate of the total mag- 
netization THcqL of the chain. In Fig. |3] we show oc/J(L) as a 
function of L for different values of A. For small L, Ji{L) oc L" 
with fl ~ 1/3 up to a length scale L* ~ O (10^) beyond which 
the effective coupling between strongly coupled clusters and 
the bath, a/7, gets larger than one and the clusters become 
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FIG. 4: a) LVp^(log(ro//i)) as a function of \og{To/h)lL^ with 
\|f = 0.31 for different system sizes for a = 0.2 and Jq = 0.025. 
b) Magnetic moment a'p{L) as a function of the system size L for 
a = 0.2 and Jo = 0.025 suggesting Tweq > and L* ~ 100. c) Ai 
as a function of L on a linear-log plot for different values of Jo and 
a = 0.2. d) Phase diagram for Iiq = 1 and Q./p = lO'* characterized 
by a single phase with nigq > 0. The dotted line represents the line of 
smeared transitions characterized by an IRFP scaling <7I9> . 



localized. Above this value /7(L) ^ Tn^qL (see inset of Fig. 
|3}, which suggests a finite magnetization Tweq before the pu- 
tative critical point is reached. This is a manifestation of the 
"frozen" clusters and lead to the concept of rounded quantum 
phase transitions in the presence of dissipation 1 14]. The typi- 
cal size of a frozen cluster turns out to be rather large L* > lO''^ 
for this range of parameters (a = 0.03 — 0.052, 7o = 0.34). 
Consequently the fraction of non-frozen samples, Ai, in Q is 
close to 1 for the system sizes that we could study numerically. 

A stronger dissipation strength a reduces L* and gives us 
the possibility to study the crossover to a regime that is dom- 
inated by frozen clusters, in particular the L-dependence of 
Al in 0, and we consider a ~ 0.2 as an example now. For 
the restricted distribution P(log(ro//i)), we obtain the same 
scenario as for smaller dissipation, as shown in Fig. |^ for the 
putative critical point Joc = 0.025. Fig. , shows (Xp{L) indi- 
cating that /J(L) o= THeqL for L > 100, which implies L* ~ 100. 
The fraction of non-frozen samples, Al, shows a clear devi- 
ation from unity already for the system sizes we study here: 
Fig. ^ shows Al as a function of L for different values of 
Jo- The data imply an exponential decay Al °^ e^^/^. The 
characteristic decay length fits well to L 7" , with (3 ~ 0.8, 
meaning a very rapid increase of L with decreasing dissipation 
strength a. By comparing L with L* for various parameters 
{a, Jo) we find that L — kL*, with k a dimensionless number 
of order one, weakly dependent on a and Jo- 

As long as L < L* the restricted distribution is not signif- 
icantly different from the full distribution of non-vanishing 
excitation energies, since the probability for a frozen sam- 
ple is small for L <C L*. Since Pl has a power law tail 
down to excitation energies exponentially small in L, the spe- 
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cific heat, susceptibility etc. in finite size systems display a 
singular low temperature behavior characterized by the dy- 
namical exponent z(a) down to very low temperatures (actu- 
ally down to Ti ~ e^"^). This intermittent singular behavior, 
X(r) - for the susceptibiHty and c{T) - r^/-'(") 

for the specific heat, persists for larger system sizes as well as 
for L — > oo, but as soon as L > L* it will eventually compete 
with the temperature dependence of the (quantum mechani- 
cally) frozen clusters - e.g. \/T for the susceptibility. Since 
the latter has a small amplitude proportional to 1/L*, classical 
temperature dependence will only set in below T* ^ 
and Griffiths-like behavior is visible (also in the infinite sys- 
tem) above T* . 

It is instructive to consider the RTIC without dissipation, 
but with a finite fraction p of zero transverse fields (i.e. p{h) = 
p5(/!) + /iq ^(1 - p)e(/i)e(/io - h)). The sites with h = Q then 
correspond to frozen clusters that have an average distance 
L* p^^ Indeed the distribution P{h/TQ) shows the same 
behavior as in Eq. with ~ e^^l^' . But, in contrast to 
the dissipative case, the restricted distribution P(log (Fo//!)) 
is naturally identical with the one for the non-diluted (p = 
0) RTIC, which shows the IRFP scaling (|9} at = 7o with 
Wtic = 0.5, different from the one we obtain here. 

The connected correlation function Cir) — 
{a]a]^^) — {&j){a]j^_i.) decays exponentially for r :s> L*, 
given that the quantum fluctuations are exponentially sup- 
pressed beyond this length scale Q, consistent with 11411 . It 
should also be noted that the connected correlation_function 
of the restricted ensemble of non-frozen samples C{r) does 
not behave critically since the number of non frozen spins 
belonging to the same cluster is bounded by 1/a in the 
restricted ensemble. Thus, the origin of the systematic broad- 
ening of the distribution P{logro/h) is here different from 
a standard IRFP and probably stems from the non-localized 
clusters with a, close to (but smaller than) one (see Eq. O. 

We have checked that the behavior of the gap distribution 
characterized by Eq. ( I7I9> depends very weakly on the ad hoc 
parameter Q./p in the range 10 — 10"*. In this range, the rela- 
tive variations of the estimated exponent \\t is of the order of 
5%, although the values of L* and (UcJoc) are more sensi- 
tive, and probably non universal. We repeated the previous 
analysis for different values of {a, Jo) (keeping ho ~ 1). In 
contrast to the pure case fT6], the entire plane is here char- 
acterized by a single phase where Tweq > 0, beyond a length 
scale L* = L*(a,J()), everywhere (except on the boundaries 
(a,0) and (0,yo)) 1 20]. One can identify a line of smeared 
transitions associated with the broadening of the restricted gap 
distribution P(log (Fq / h)), according to (|9j: this is depicted in 
Fig. We find that the associated exponent \\t vary weakly 
along this line, its relative variation being less than 10%. 

To conclude our strong disorder RG study of the RTIC cou- 
pled to a dissipative Ohmic bath revealed that non-frozen sam- 
ples display an IRFP scaling of the distribution of excitation 
energies. With this we computed a continuously varying ex- 
ponent z(a) that determines an intermittent singular tempera- 
ture dependence of thermodynamic quantities above a temper- 



ature T* ~ L*^'"("^. L* is a characteristic length scale above 
which the ground state displays a non-vanishing magnetiza- 
tion, as predicted by the smeared transition scenario 1 14|, and 
that we determined to increase exponentially with the inverse 
strength of the dissipative coupling. This implies that numer- 
ical studies can hardly track the asymptotic behavior JitIi and 
that experiments at very low but non-vanishing temperatures 
might still show indications for quantum Griffiths behavior 
llOl I12I I13I1 . In higher dimensions we expect a similar sce- 
nario as the one discussed here and it would be interesting 
to extend our study to Heisenberg antiferromagnets and XY 
systems. 
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